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Coulomb and mass difference corrections to fow-energy pion-nucleon scattering are calculated 
and compared with previous work including potential models, dispersion relation methods and 
chiral perturbation theory calculations. Particular attention is paid to their role in testing isospin 
breaking. 



I. INTRODUCTION 

Electromagnetic corrections are a key element in understanding low energy pion nucleon scattering and, in 
turn, the fundamental elements of chiral symmetry breaking and the pion-nucleon coupling constant. They 
are also crucial in the important role of low-energy pion-nucleon scattering as a testing ground for isospin 
breaking. 

Tests of isospin breaking have been re portedQ, H i using the extracted charged pion scattering ampli- 
tudes to predict charge exchange amplitudes. A breaking of the order of 8% in the amplitude was found. 
Piekarewicz 4] has suggested a possible explanation for the isospin breaking based on the effect of quark 
mass differences on the charged and neutral pion coupling constants. 

The breaking mentioned above is beyond the expected corrections due to electromagnetic interactions 
and mass differences. Corrections for these effects have been reported using a potential model for their 
calculation^ ^ based on pioneering work by Oades and RacheQ. In addition, an analysis was madeQ 
(hereafter referred to as GAK) in which the potential method corrections were incorporated into the fit itself. 
These last corrections have not been previously reported and it is a primary purpose of this paper to give 
them. 

For unbroken isospin one can write the charged pion scattering amplitudes in terms of the isospin ampli- 
tudes. 



a) 

In the potential treatment each of the pure isospin amplitudes is assumed to satisfy a wave equation 
with potentials corresponding to isospin [V^ and V'^). Since the Tr+p amplitude consists of a single isospin 
amplitude, only the Coulomb correction needs to be taken into account. The method for this case is covered 
in Section ITU 

The last two amplitudes, being a mixture of isospin amplitudes, pose the additional problem of the 
neutron-proton and tt^-tt*^ mass differences. By using a linear combination of the amplitudes, the solution 
is expressed as a pair of coupled equations corresponding to the charge states. 
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T.Vc + (^Ivi + i^i) V, + ^ (^^ - ^0 = ec^c (2) 

ToV'o + ^ (^^ - V^^) iPc + (^V^- + "^Vi^ iPo = eoVo (3) 

where ipc and ec are the charged wave function and kinetic energy and tjjo and eo are the corresponding neutral 
pion quantities. Tc and To are the kinetic energy operators. The corrections are then included by adding 
a Coulomb potential to the equation describing the tt^ only. The masses entering into the kinetic energy 
operators, T, and the kinetic energies, e, are also replaced by expressions with the true masses corresponding 
to the charged or neutral pions and nucleons. It is assumed that the potentials remain unchanged during 
this process. This assumption has been questioned by RusetskyQ. 

Although the general method used in each case is very similar, each group has used different forms for the 
potentials, fit different data sets, used different relativistic prescriptions for the pion energies and reduced 
masses (energies) and, in some cases, different wave equations. Thus, variations are to be expected among 
the resulting corrections. 

In addition to the potential method, a dispersion approach has been used Sauter^ developed 

the basic method for s-waves only without considering mass differences. The model was extended in several 
works fT^. This approach has become the most common method for making the corrections and they are 
often referred to as the NORDITAIH ini results. The advantages and disadvantages of the potential and 
dispersion relations method are discussed in Ref. I^. 

II. 7r+ PROTON SCATTERING 

In this case the corrections in a potential model are straightforward. The wave equation can be solved for 
each partial wave to obtain the phase shifts. While this full solution for the amplitude includes everything, 
the partial wave expansion diverges, as does the one for a pure Coulomb potential. Since the sum of the 
Coulomb series alone is known, the non-spin-flip and spin-flip amplitudes can be written as 

f{9) - ,U9) + -1- + l){e^'''+ - e^-^ + ^(e^*- - e^'^')]Pe{co8 9) (4) 

5W = ^Et^'^''" - e'*-]P/(cos^) (5) 

where ag is the Coulomb phase shift and /c(£) is the Coulomb amplitude. The notation £zL means j = £zt^. 
Here (and in GAKQ) (unlike Gashi et al.[3) we neglect the spin dependent part of the electromagnetic 
interaction. Eq. ^is often written as 

oo 

fie) = M9) + ^ E [(^ + - 1) + ^(^"'""^ - (6) 

where the nuclear phase shift is defined as 

- 6,± - a,. (7) 
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FIG. 1: Comparison among the different method for n'^p elastic scattering. The difference between the point charge 
(pc) and full calculation of Gashi et al.Q is very small except for the s wave. The dotted curve represents a fit to 
the NORDITA resultsilflj to be used later and given in the appendix. 



This phase shift is then compared with the hadronic phase shift, d'^_^_, i.e. the one that would exist if there 
were no Coulomb interaction. A potential model is used to calculate these two quantities in the present case. 
In the limit that the interactions are very weak, the hadronic and nuclear phase shifts become equal so that 
the difference may be expected to be small. The correction is defined as the difference between them. 



Ci± = (5"± - (8) 

Figure ^ shows a comparison of the results the different determinations. It is seen that the potential 
methods are in reasonable agreement. The point Coulomb result of Gashi et al.Q included no higher order 
electromagnetic effects whereas the full calculation included a finite charge distribution as well. Since the 
corrections from GAK included a finite charge density but no no higher order electromagnetic effects, there is 
no exact comparison possible. The dispersion relation approach is seen to give only slightly different results. 

While one can judge to some extent the size and degree of agreement from Fig. ^ it is useful to see the 
effect on the cross section itself. Figure shows such a comparison at a kinetic energy of 45 MeV. The 
hadronic phase shifts are taken for Ref. 1.] for all of the cases. It is seen that the correction is small but 
not negligible, being of the order of two standard deviations of accurate data. The differences between the 
corrections would seem to be unimportant at the present level of the data. 
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FIG. 2: Effect of the corrections for tt p scattering in a typical case. It is seen that the corrections are small but not 
negligible. The three methods compared give very nearly the same result. The data are from Brack et al. [Tsf . 



III. TT PROTON SCATTERING AND CHARGE EXCHANGE 



In this case the expression of the corrections is more complicated since the potential system consists of 
a pair of coupled differential equations. Not only is there a correction due to the finite charge in the Tr^p 
channel but the neutron-proton and 7r° — tt^ mass differences need to be taken into account. 



A. Method 



To express the corrections in the coupled channel case we follow the formulation of Gashi et al. 

i- The 

solution of the coupled channel system^, 0| produces the amplitude for tt^ proton scattering, charge 
exchange and tt" neutron scattering. The result is expressed as a 2 by 2 symmetric S matrix with the diagonal 
elements representing 7r~p and Tr^n elastic scattering and the off-diagonal elements charge exchange. The 
scattering matrix is expressed in terms of a 2 by 2 real K matrix. 



l + iK .S-l 

These equations hold for each partial wave where we have suppressed their indices. The Sn matrix element 
has been multiplied by e"^'"'* to remove the pure Coulomb effect before the calculation of the K matrix. 
Since the K-matrix is real and symmetric it can be transformed to a diagonal form with an orthogonal 
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FIG. 3: The solid curvegives the result with both Coulomb and mass difference corrections from GAK. The dashed 
line is from Gashi et al.jg], the solid dots are from Tromborg et al. lOj and the solid triangles are from Zimmermann 7]. 
The dotted hnes represent a fit to the NORDITA 10] results used later. 



If isospin is a good quantum number then cos — > y | = cos 0o a-nd <5i and 82 are the isospin i and | 
phase shifts. The form expressed in Eq. 1101 is valid independent of isospin conservation and the differences of 
the resultant phase shifts 5i , 62 as well as the value of (j) needed to diagonalize the system, from the unbroken 
isospin values are used to quantify the breaking. 

The corrections are expressed as 



22 22 



(11) 
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FIG. 4: Comparison with the correction for the rotation angle Atj). The meaning of the symbols is the same as in 
FigH 



for each partial wave, where the barred quantities are those obtained from the model by solving for Si , 62 and 
(j) from the equations above with Coulomb and mass differences present. The quantities with super script 
"h" (the hadronic values) are those obtained from the solution of the system with all pion masses equal to 
the charged pion mass, the neutron mass equal to the proton mass and no Coulomb interaction. 



The hope is that these corrections will be nearly model independent and one can fit the data by choosing 
S'l and 62, apply the corrections to obtain the full K and S matrices, add in the Coulomb amplitude where 
appropriate and compare with data to calculate a value of x^. A search on the values of the isospin pure 
phase shifts, 61 and 62 can then be made. 

We now address the problem of calculating the corrections using the model. Since Ki and K2 are the 
eigenvalues of the K-matrix we can solve directly for them with 

JTii + K22 ± ViKii - K22Y + ^K'i^ 
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FIG. 5: The short dashed hne contains the Coulomb correction only, the dash-dot line is with the mass difference 
correction only and the solid line contains both corrections. The results of Zimmermann are given by the solid 
triangles (full calculation) and solid squares (Coulomb correction only). 



We may now identify the plus and minus signs in this expression with the isospin i and | states. 

Multiplying Eq. 1101 from the left by the transpose of the rotation matrix we find four relations for the 
tan(/<, all of which are equivalent. Two of them are: 

- Ki Ki2 , , 

tan0= ; tan0=— — . (13) 

Ki2 K2-K11 



B. Difficulty with the Pj. wave 

While in most cases there is no difficulty in choosing the proper sign in Eq. 1121 to make the correct 
association of and 82 with the isospin phase shifts, in the case of the Pi scattering a problem arises. In 

this case around 60 MeV the two hadronic phase shifts (P31 and Pn) cross. For the two eigenvalues, Ki 
and K2 to be equal the radical in Eq. II 21 must vanish. This requires that the conditions 
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FIG. 6: Individual contributions to A(j>. The convention for the lines and points is the same as in Fig. |3 



be satisfied simultaneously. For unbroken isospin tliis indeed happens. In this case isospin requires 

A'i2 = V2{K22 - Kii) (15) 

so that 



A 



(16) 



In order for A to have a continuous derivative the ± association must be changed at the zero of the radical. 
If one chooses the sign such that the radical has the actual sign of /Cn — A'22 then 



A+ = 2Kii - K22 = tan (5 2 ; A_ = 2K22 - Kn = tan (5 3 



(17) 



For broken isospin, in general, both conditions in Eq. ^^virill NOT be satisfied at once so that the two 
eigenvalues can never be equal. To understand what happens in this case it is useful to look at the direct 
solution of Eq. ^| 

In the general case we can expand Eq. ^| to find that the condition for the off-diagonal elements to be 
zero is 



2(j) = tan 



_1 2^12 

^^22 - Kn 



(18) 



9 



160 h I I I I I H 




40 F ^ 

10 20 30 40 50 60 



FIG. 7: Com pari son of methods of corrections for charge exchange for forward angles at low energies with the data 
of Frlez et al.|l6l|. The dotted curve represents no correction, the solid curve gives the correction from GAKQ, the 
chain-dash curve is from Gashi et al.0 and the dashed curve represents the corrections from NORDITA[io) 



so that 20 changes by tt around the point K22 = Kn. Consider a simple example in which the isospin 
relation is slightly broken in a simple way such that 

K12 ^ V2{K22 - Kii) + e (19) 

then 

2(t> = tan-i {2V2 + — ) (20) 

\ A22-Aiiy 

so that 20 increases by tt as the difference K22 — Kn changes sign. We have assumed that K22 > Kn for 
energies lower than the crossing point. If the opposite is true, a very similar argument gives the same effect 
with 2(f> decreasing by tt. 

The values of Ki and K2 are given by 

Ki = iii'ii(l + cos2(/)) + ix22(l-cos20)-Xi2sin20'^^*''tan(S5 (21) 



K2 ^ -is:ii(l-cos20) + -X22(l+cos20) + Xi2sin20'^^*''tan(S5 (22) 

Indeed, as 20 increases by tt, Ki ^ K2 since the sine and cosine change sign under these conditions. 

Thus, the correction is non-perturbative. The change in 20 by tt occurs at the point where the hadronic 
phase shifts cross and hence depends strongly on the hadronic phase shifts used to define the corrections. In 
this case the general method fails except for corrections far from the crossing point. 
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One can choose to change the association of the plus and minus signs at the zero of K12 (that is what is 
done in the present work) or when Kn = K22 but either way there wiU be discontinuities in the corrections 
extracted at the crossing point. The actual amplitudes remain continuous, it is only this representation of 
the breaking which shows an anomalous behavior but the compensation only occurs for the hadronic phase 
shifts used in the calculation of the correction. If one uses these corrections with phase shifts other than 
the ones for which the corrections were derived, discontinuities in the resulting amplitudes will result. Of 
course, this is a small correction in a small partial wave in this case so in a practical sense the problem may 
not be very serious. 

The difficulty comes partly from the fact that, even for good isospin, for Ki = K2 any value of (j) is valid for 
the transformation so is undetermined at that point. For unbroken isospin there is no problem, however, 
since is determined arbitrarily close to the crossing point on either side. 



C. Comparison with previous models 

An obvious isospin breaking effect due to the mass differences leads to a positive Q-value for n^p — )■ n'^n 
and a negative Q-value for n^n —> Tr~p. This fact destroys the symmetry of the amplitudes derived from 
the symmetric S-matrix, due to a factor kf/ki originating from the incoming and outgoing fluxes which is 
obtained in the charge exchange cross section obtained from the procedure above. This factor is included in 
the calculations of charge exchange cross sections to be shown in this section. 

Figures 131 and ^ show a comparison of the corrections obtained from GAK 1] with previous work. It is 
seen that there is qualitative agreement among all calculations. There is rather good agreement between 
the present work and Zimmermann at the higher energies where he calculated. For the Sn and P33 partial 
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FIG. 9; Comparison of the effects of the corrections with the charge exchange data of Isenhower et al.[i7 
MeV. Other comparisons at 10.6 and 20.6 MeV are similar. 



at 39.4 



waves there is good agreement with the work of Gashi et al.|3 as for the case of A0 for the Ps channel 
where, in fact, all calculations agree. 

Figures El and [HI give the breakdown of the contributions to the corrections into Coulomb and mass 
differences for GAK. It is seen that the most important contribution is from the mass differences, especially 
at low energies. Since Zimmermann gave a separation of the two effects for the P33 phase shift and the 
A(/) for the P| channel a comparison can be made and the agreement is good. The NORDITA^^I group 
also gave a separation into the Coulomb and mass difference contributions. The mass difference correction 
dominates at low energy (see Ref. Figures 2, 3 and 4). Hence, while these corrections are often called 
electromagnetic, it is the correction arising from mass differences which is the most important. 

While the expression of the correction in this form is very useful, it is difficult to appreciate the effect 
on the cross sections of the differences among the determinations as they may well be compensating to give 
similar corrections to physical observables. Hence, we consider the effects of these corrections on some charge 
exchange cross sections and the 7r~p elastic amplitude. In each case the "no correction" calculation is made 
from the hadronic phase shifts of GAK^I and the corrections are made relative to them. 

Figure [7| shows a comparison of the correction with the charge exchange data of Frlez et al.^^ at 27.5 
MeV for the three groups who calculated in this energy range. It is seen that the resulting corrected cross 
sections are all very nearly the same and that the correction takes the result away from the data whereas the 
uncorrected calculation (dotted line) lies much closer to agreement. This last effect occurs only at forward 
angles. 

Figure IHl gives the breakdown of the contributions. It is seen that the Coulomb correction actually moves 
the result slightly closer to the data while the mass difference correction takes it away. 

Figure m shows the effect over the full angular range compared with the data of Isenhower et al.0|. It is 
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FIG. 10: Comparison of the effect of the corrections with the data of Bagheri at al.[l^ at 62.2 MeV. Other comparisons 
with data from the same group at 45 and 92 MeV are similar. 



seen that for larger angles the corrections move the prediction closer to the data. 

At higher energies the effect of the corrections is somewhat smaller. A comparison with the data of Bagheri 
et al.[l^ is shown in Fig. ^| The experimental results were given as coefficients of Legendre polynomials 
in this case and the dots shown were calculated from these numbers. The effect of the corrections is seen to 
be small at this energy and the two potential models agree well. 

In general, the difference between prediction and data can be expressed as a difference in normalization 
of around 15%. 



The forward charge exchange cross section is a sensitive measure of isospin breaking since it shows a deep 
minimum near 45 MeV. While most checks on isospin depend on obtaining an absolute cross section, the 
position of this minimum does not; only a knowledge of the beam momentum is needed. 

The left part of Fig. ^2 shows the effect of the corrections of the three groups compared with the data of 
Fitzgerald et al.[3- We see that the corrections are all very nearly the same and move the predicted cross 
section away from the data. While the uncorrected prediction cannot be said to give a good fit to the data, 
it is much better than the prediction after correction for the mass differences. The Coulomb potential has 
very little effect on the position of the minimum. 

The right part of Figure ITTI shows the prediction of the corrections for the amplitude for ir^p n^p. It 
is seen that the mass correction is significantly different from the pure Coulomb effect. One can also note a 
large difference among the methods. 



13 




FIG. 11: Left: Zero degree charge exchange with corrections from GAkQ, Gashi et al. andNORDITA[l3. Right: 
Corrections to the real part of /,r-p->7r-p GAKQ (dash: Coulomb only, solid: full correction), Gashi et al.0 
(dash-dot) and NORDITAtlfl| (dotted). 



D. Comparison with Chiral Perturbation Theory 



While the results of the previous section are rather consistent, at least for the charge exchange, they may 
all be wrong. Thus, it is very desirable to have a comparison with an independent method calculating from 
an approximate QCD viewpoint. 

A calculation of isospin breaking in Chiral Perturbation Theory (ChPT) has been made by Fettes and 
MeiBner 20). They calculated 6 ratios which are measures of isospin breaking. Most of the ratios (3-6) involve 
quantities difficult to measure. Ratio number one involves isoscalar quantities which vanish at threshold so 
that the ratio will depend sensitively on the phase shift fit used to compute them. 

Their second ratio, 

' = f TT TTTf ^ ^ 

involves the same expression as is used to test isospin in the charge exchange reaction so is very suitable for 
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comparison. As stated by the authors, a direct comparison can not be made between their work and data 
but one can calculate the ratio with the hadronic model under the same conditions. 

If one were to attempt to compare with actual charge exchange amplitudes obtained directly from the 
data they would include the flux factor ^JkJJki. Since the ChPT calculation is made at the matrix element 
level, this factor should not be in the comparison. We show both cases to demonstrate the importance of 
the factor. 

Since they calculated corrections to the strong process their combination 

goes to zero in the limit of no strong interaction and we must remove the pure Coulomb amplitude. We might 
subtract it partial wave by partial wave or simply use the nuclear phase shifts. It is the latter prescription 
which is shown in the comparison graphs. There is, in fact, very little difference in these two ways of removing 
the Coulomb effect. 



Figure IT^ shows the comparison with and without the flux factor. We see that, for the S and Pi waves 
the pure Coulomb correction agrees rather well with the ChPT (full) calculations. The Pi ratio is large and 
shows a singularity around 60 MeV due to crossing of the hadronic phase shifts leading to the vanishing of 
the charge exchange amplitude in this partial wave. The position of this crossing depends on the detailed 
fit of the hadronic phase shifts and hence cannot be expected to be the same in the different cases. For the 
Pa wave, even the ratio calculated with Coulomb corrections alone is more than a factor of 2 larger than 
the ChPT result, perhaps due to the fact that the amplitude is larger and the perturbation series has not 
converged. 

Without the flux factor, only the Pa partial wave shows a significant discrepancy between the models and 
the ChPT calculations. The two potential models and the dispersion relation approach show about the same 
factor of 3-5 difference in this case. Since this is the channel in which the Delta resonance occurs, one can 
question whether or not the explicit degrees of freedom of the Delta can make an important difference. 



IV. CONCLUSIONS 

We have seen that the hadronicly calculated corrections of GAK[l|are in rough agreement with previous 
determinations. In most cases, the values of GAK and Gashi et ai.0, 1^ are in reasonable agreement. All 
of these methods give a significant (and very similar) correction to the charge exchange cross section at low 
energies due to the mass differences. 

Corrections calculated from ChPT are significantly smaller than those of the hadronic methods in the Ps 
partial wave, especially when the mass correction is included. This may indicate that there is a general error 
in the hadronic corrections or that the Delta degree of freedom needs to be treated differently in the ChPT 
calculation. 

Since the mass correction in the form currently used has been questioned it is possible that the assumption 
about the isospin potentials remaining unchanged is incorrect. 

If one assumes that the hadronic corrections are too large then a closer agreement can be achieved in the 
case of forward angle charge exchange at low energy and the position of the forward minimum. However, 
this conjecture does not correct the discrepancy at larger angles. So even this rather drastic possibility does 
not solve the problem. 

We thank G. C. Oades and Ulf Meissner for comments on the manuscript. 

This work was supported by the National Science Foundation under contract PHY-0099729. 
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FIG. 12; The left figure sliows tlie ratio R2 calculated with GAK corrections for Coulomb only (dash-dot) and 
Coulomb plus mass differences (solid) compared with that from Fettes and Meifiner.SOj (solid squares). The dashed 
line gives the result with the corrections of Gashi et al. 6] and the dotted line gives that of NORDITA 10] . The flux 
factor y^kf /ki is included in the hadronic calculations. The right flgure shows the same information except that the 
flux factor \Jkj jki is not included in the hadronic calculations. 



APPENDIX A: FITS 

Following are a set of fits to the phase shifts and corrections which are given as an aide. The form of the 
fits has no physical significance and should not be used outside the range 5-100 MeV. In what follows, T is 
the laboratory kinetic energy in MeV and 8 is in degrees. 



1. GAK Hadronic Phase Shifts 



The fit to the GAK hadronic phase shifts for the S waves is in the form 

1 + 5T + cT2 ■ 



(Al) 
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For Si a = 1.0319; b = 2.865 x 10"^; c = -8.0 x 10"^ (A2) 

For S3 -0.5197; b = -5.37 x lO"'"^; c = 1.939 x 10"^ (A3) 

For the p-waves 

0.001806Ti(l - 0.00414r) , 0.01084Ti(H- 0.0085T)(1 + .06e-° ''"7') 

01.3 = ; 033 = A4) 

1 + O.OOIT ' 1 - 0.000436r ^ ' 

0.00636ri(l - 0.006190r- 0.0304r2) ^ 0.00196ri(l + 0.00228T) , 

Oil = ; ^3^ — A5) 

1 + 0.0131T ' l + O.OOOllT ^ ^ 

2. TT+p Corrections 

The fits to the 7r+p corrections by GAK are 

Co+ = 0.1 + 0.000049r; Ci_ = -0.009 + 0.00045T; Ci+ = -0.0185 - 0.000397ri-5 (A6) 
The fits to the 7r+p corrections by Gashi et al.0 are 

Co+ = 0.075 + 0.0005r; Ci_ = 0.005 + 0.0002T; Ci+ = -0.033 - O.OOOOIQST^ ^ (A7) 
For NORDITA 

Co+ = 0.095; Ci_ = 0.00045T; Ci+ = -0.035 - 0.000067r2 (A8) 

3. Coupled Channel Corrections 

The fit to the GAK corrections for the s wave is 

C,= 0.753T-ni-0.0114T); C3 ^ - "'^^^^"^ ' "'""^^^^ (A9) 
^ ^' l.- 0.00299r ^ ^ 

A(/)= ^+0.0011T- 0.026 (AlO) 

For the partial wave 

Ci = -0.035- 0.01745 tan-^ f^:^) (All) 

1. + .0008T2 Vr-60.1/ ^ ' 
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01137T'' '^^^ / 0155 \ 

- - (1 + 0.137TO.^^ + 0.000039T^) ~ (7-60! ) ^^^^^ 

For the Pa partial wave 

Ci = -0.001-0.0039r° ''-0.0038 f — j ; C3 = 0.1+0.0048T-0.23 f — j ; A0 = — ^+0.37 (A14) 



The fit to the Gashi et aL0 corrections used in this work is 



Ci-H^- 0.056; C3--^- 0.04 (A15) 



For p| Ci = ~0.0005T + 0.006; C3 = 0.0036T + 0.13 (A16) 
For pi Ci = -0.045 + 0.000016Ti-^; C3 = -0.05 + 0.00025T + 0.046e"° °^^ (A17) 



A'^^ = ^ + 0-0013r- 0.149; A0p, = - 1.9; A</.p3 - + 0.4; (A18) 

The fit to the N0RDITA[i3 corrections used in this work is 

Fo^^ Ci = i^^j - 0.095; ^3 = -^^ + 0.082 (A19) 



For p3 Ci = -0.00065r; C3 = 0.008T + 0.06 (A20) 



= ^ + 0.0008T + 0.01; A0p^=-^ + O.4 (A21) 
The NORDITA group did not give corrections for the Pi partial wave. 
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